The Ekman-Hartmann layer in MHD Taylor-Couette flow 
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We study magnetic effects induced by rigidly rotating plates enclosing a cylindrical MHD Taylor- 
Couette flow at the finite aspect ratio H/D = 10. The fluid confined between the cylinders is 
assumed to be liquid metal characterized by small magnetic Prandtl number, the cylinders are 
perfectly conducting, an axial magnetic field is imposed Ha « 10, the rotation rates correspond 
to Re of order 10 2 — 10 3 . We show that the end-plates introduce, besides the well known Ekman 
circulation, similar magnetic effects which arise for infinite, rotating plates, horizontally unbounded 
by any walls. In particular there exists the Hartmann current which penetrates the fluid, turns into 
the radial direction and together with the applied magnetic field gives rise to a force. Consequently 
the flow can be compared with a Taylor-Dean flow driven by an azimuthal pressure gradient. We 
analyze stability of such flows and show that the currents induced by the plates can give rise to 
instability for the considered parameters. When designing an MHD Taylor-Couette experiment, a 
special care must be taken concerning the vertical magnetic boundaries so they do not significantly 
alter the rotational profile. 

PACS numbers: 47.20.-k, 52.30.Cv, 47.15.Cb, 52.72. +v 



I. INTRODUCTION 

Motion of a fluid confined between two concentric, ro- 
tating cylinders is a classical problem in hydrodynam- 
ics and, if the fluid is conducting and an external mag- 
netic field is applied, magnetohydrodynamics (MHD). 
The flow of this type, usually referred to as the Taylor- 
Couette flow, has been first studied by Couette JM and 
later was subject of a seminal work by Taylor [2J, who 
experimentally confirmed theoretical results of a linear 
stability analysis. In the field of MHD, an important 
work was done by Velikhov [|[ who has shown that for 
the conducting fluid a weak magnetic field can play a 
destabilizing role and can lead to an instability which 
today is called magnetorotational instability (MRI [4|]). 

When studying the Taylor-Couette system it is com- 
mon to assume some simplifications, the small gap ap- 
proximation or large aspect ratio. In the former it is as- 
sumed that the gap between the cylinders D = i? ou t — -Ri n 
is small compared to the radii, i.e., D/R out <C 1, this al- 
lows the neglect of terms of order 1/R, R being distance 
from the center of rotation. When considering the large 
aspect ratio, one assumes that the height of the cylinders 
H is much larger than the gap width T = H/D ^> 1, 
which guarantees that a secondary flow due to the plates 
bounding the cylinders is insignificant and does not dis- 
turb the rotational profile of the fluid. 

On the other hand, there is also plenty of work done for 
small aspect ratio r « 1, where the rigidly rotating end- 
plates play crucial role and simply introduce a new class 
of problems. When T becomes an important parameter 
it is possible to observe a wide family of different states 
(including non-axisymmetric ones or peculiar asymmetric 
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patterns - anomalous modes) for the same parameters, so 
that the observed results depend on their path through 
the parameters space from an initial state. Therefore this 
system is an excellent subject to the bifurcation theory 

In the present work we focus on the case of wide gap 
i?i n /i? out — 1/2 and T = 10 which is an intermediate as- 
pect ratio, between very short and long containers, yet in 
purely hydrodynamical contest the influence of the ver- 
tical boundaries is small, at least for Reynolds numbers 
of order O(10 2 — 10 3 ). However, if the rotation rates are 
large enough, so that the corresponding Reynolds num- 
ber is O(10 5 ) and larger, the plates can easily dominate 
the flow in the entire container. This is due to the Taylor- 
Proudman theorem, from which follows that in rapidly 
rotating systems the flow tends to align itself along the 
axis of rotation. For such rotations it is necessary that T 
would have to be several thousand in order to obtain the 
rotational profile which is not profoundly altered by the 
end-plates fill ] . 

Results of a recent MRI experiment PROMISE [H, 
EH, [HI i as well as nonlinear simulations [HI, EH indicate 
that for a flow with relatively small Reynolds number 
« 10 3 , and parameters resembling essentially MHD sta- 
ble flow in the limit of infinitely long cylinders, there ex- 
ist unexpected time-dependent fluctuation of the velocity 
field. These disturbances arise as an effect of the verti- 
cal boundary conditions, moreover the simulations show 
that they are much stronger if the end-plates bounding 
the cylinders are assumed to be perfectly conducting. 

The plates induce a well known hydrodynamical ef- 
fect - the Ekman circulation, which is a result of unbal- 
anced pressure gradients in vicinity of the vertical no-slip 
boundary conditions. There the Ekman layer develops 
in which the fluid velocity from the bulk of the container 
must match the velocity imposed by the end-plates. 

It seems that for MHD Taylor-Couette flow, magnetic 
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effects, unfike the classical hydrodynamical Ekman layer, 
induced by the plates have been overlooked. In this pa- 
per we argue that the rigidly rotating plates together 
with an imposed axial magnetic field give rise to a similar 
layer which develops for an infinite, rotating plate serv- 
ing as a boundary for the conducting fluid. One of the 
most important features of such flow is the existence of 
the Hartmann current (absent in the conventional Hart- 
mann problem, which leaves the boundary layer and 
then interacts with the magnetic field. In particular, this 
becomes important for conducting plates which was the 
case for the PROMISE experiment, since one of the end- 
plates was made from copper. 

We discuss properties of Ekman-Hartmann layers for 
infinite, rotating plates and relate it to the end-plates 
enclosing the cylinders in a Taylor-Couette setup. It is 
shown that for considered radial boundary conditions, 
the induced current turn eventually in the radial direction 
and acting in concert with the imposed axial magnetic 
field gives rise to a body force. 

We demonstrate that magnetic effects induced by the 
end-plates enclosing the cylinders can profoundly al- 
ter flow properties. In particular the rotational pro- 
file can become significantly different from the expected 
parabolic Couette solution. Moreover, if the Hartmann 
current is strong enough, it is likely that the local 
Rayleigh criterion for stability will be violated and the 
flow becomes centrifugally unstable. In an MRI exper- 
iment it is crucial to rule out such instabilities and a 
special care concerning the vertical boundary conditions 
is needed in order to obtain the desired rotational profile. 



II. PROBLEM FORMULATION 

We consider two concentric cylinders with radii 
-Rinj-Rout embedded in an external axial magnetic field. 
They rotate with angular velocities Oi n ,J7 out , the ra- 
dius ratio is fj = Ri n /Rout, the rotation ratio ft — 
^out/^in- Cylindrical coordinates (R, <j>, z) with unit vec- 
tors e^,e^,e z are used. If the cylinders are unbounded, 
i.e., infinitely long or periodic, the rotational profile is 

r> (i?) = a+A j (1) 

with 

a = ft inT -4 b=^—^Rf n Q in , (2) 
1 — ?y z 1 — r] z 

and ur = u z = everywhere. The flow is hydrodynam- 
ically stable if the Rayleigh criterion d(R 2 Q) 2 /dR > 
is fulfilled, i.e., for ft > if . Consequently, for the con- 
sidered radius ratio fj = 1/2, the flow is always stable if 
ft, > 0.25. Here we consider only cases when ft > 0.25 so 
that hydrodynamical instabilities are ruled out. 

Let us introduce the Reynolds number Re, which mea- 
sures the rotation rates, and the Hartmann number Ha, 



which measures the strength of the externally applied 
magnetic field B = B e z , 

I D 2 MinRinD 

Ha = B a J , Re= 3 

where p is the density, v the kinematic viscosity, rj is 
the magnetic diffusivity, /io is the magnetic permeability. 
The fluid confined between the cylinders is assumed to be 
incompressible and it can be characterized by the mag- 
netic Prandtl number Pm = vjr\. For laboratory liquid 
metals, like gallium, Pm is very small - of order 10~^ 5 " 6 \ 
therefore we concentrate on effects arising only when Pm 
is small. 



A. The Equations 

Using D as the unit of length, v/D as the unit of ve- 
locity, D 2 jv as the unit of time, Bq as the unit of the 
axial magnetic field and assuming B = Bo + b we can 
write non-dimensional MHD equations for the problem 
of our interest, i.e. 

d t u + (u ■ V)u = -Vp + V 2 u 

Ha 2 

+ TT~ [(rotb) x b+ (rotb) x B Q /B Q ] , (4a) 
Pm 

d t b = 7^-V 2 b + rot(u x b) + rot(u x B /B ) 7 (4b) 
Pm 

with divu = divb = 0, where u and b are the velocity 
and the perturbed magnetic field, p is the pressure. 

For the velocity we apply no-slip boundary conditions 
at the cylinders and at the end-plates as well. We assume 
that both the plates rotate rigidly with angular velocity 
^ond, which can be set to any value so that the plates 
can rotate independently of the cylinders. 

Boundary conditions for the magnetic field are deter- 
mined by magnetic properties of the cylinders and the 
plates. Here we consider only perfectly conducting ra- 
dial boundaries, so that the transverse currents and per- 
pendicular component of the magnetic field vanish, hence 
R-^+dRbj, = at R = R in /D, R = R out /D. We chose 
such boundaries since in the PROMISE experiment the 
cylinders were made of copper. The reason for choos- 
ing copper is that the critical Re and Ha numbers for 
the onset of the MRI are smaller by almost a factor of 2 
with perfectly conducting boundaries than with insulat- 
ing boundaries (l8| . 

For the end-plates, similarly like for the walls, the elec- 
tric field must be continuous and b z — 0, then b^ = ed z b,p 
at z = and 6^ = —td z b^, at z = T where e characterizes 
a thin layer of relative conductance of the fluid and the 
plates 0, |2(j. When e — > we obtain conditions cor- 
responding to insulating end-plates, i.e., b$ = d z j$ = 
at z = 0, z — r, jt/, being the azimuthal current. For 
e — > oo we have the case describing the perfectly con- 
ducting plates, d z b^ — j$ = at z = 0, z — F, We note 



3 



that this thin-wall approximation is valid only when the 
magnetic field varies linearly within the plates and it does 
not necessarily resembles situation in a real experiment. 

B. The small Pm limit 

For laboratory liquids the conductivity a is small, so 
that the magnetic diffusivity r\ = l/^o " is very large 
(compared to the viscosity) and the corresponding mag- 
netic Prandtl number Pm is small. Consequently the 
time scale for magnetic diffusion is much shorter than 
other time scales. Therefore we consider the limit r\ — > 
oo, however it must be supposed that Ha tends to a fi- 
nite value. The perturbations b of the externally applied 
field induced by the motion of the fluid are Pm times 
smaller than Bo, although theirs effect on the Lorentz 
force can not be neglected since Ha 2 /Pm[(rot b) x Bo/ -Bo] 
is already of order Ha. Nevertheless the interactions 
(rot b) x b are vanishingly small. 

Similarly in the induction equation we may apply a 
quasi-static approximation, so that the electromagnetic 
field proceeds along a sequence of steady-state solutions 
of the Maxwell equations to conditions described by u, 
and therefore b in each moment adjusts instantaneously 
to the velocity u. Hence, in the small Prandtl limit 
Pm — > 0, the system ((4]) can be written as 

d t u + (u • V)u = -Vp + V 2 u+ 

Ha 2 (V x b) x Bo/So, (5a) 
V 2 b = -V x (u x Bo/So) (5b) 

with divu = divb = 0, [HI, HI- The equations © to- 
gether with the discussed boundary conditions are solved 
with finite difference method using stream function- 
vorticity formulation in the (R, z) plane. In this work 
we assume that the flow is axisymmetric. For more de- 
tails on the numerical procedure see [l(| HH • 

III. THE EKMAN-HARTMANN LAYER IN THE 
MHD TAYLOR-COUETTE FLOW 

At an interface between an incompressible fluid with 
low viscosity and a rapidly rotating rigid surface devel- 
ops an Ekman layer with thickness oc v/£l, where 
n is rate of uniform rotation. Similarly for a flow of con- 
ducting, incompressible fluid in vicinity of a rigid non- 
rotating boundary, and under the influence of an external 
magnetic field perpendicular to the surface there exists a 
Hartman layer with thickness c?h oc Ha -1 . When these 
two effects are combined, the Ekman-Hartmann layer de- 
velops 23]. It can be viewed either as a modification of 
the Ekman layer by introducing the conducting fluid and 
imposing the external magnetic field or as a modification 
of the Hartmann layer by adding the uniform rotation of 
the bounding surface. The resulting layer (in its steady 
form) assures a proper transition for the velocity and the 



magnetic field from values inside the bulk of the fluid to 
the applied boundary conditions. 

The linear analysis of the Ekman-Hartmann layer in its 
idealized case was presented by Gilman and Benton [24[ • 
They have considered an infinite, insulating plate rotat- 
ing with Opiate at z = 0, and a conducting fluid filling 
the space z > 0, the fluid far from the plate rotates with 
^fluid = Opiate (1+e), e C 1. The most important conclu- 
sion of this work was that in addition to the well known 
Ekman suction/blowing of mass flux there also exists an 
electric Hartmann current which has the same direction 
(or opposite when the external B z is negative) as the ve- 
locity of the Ekman blowing (that if fluid is blown away 
or sucked towards the boundaries depends only on sign 
of e). This current, which is arising due to the vertical 
shears, leaves the Ekman-Hartmann layer and potentially 
influences the flow far away from the boundary. 

For the magnetized Taylor-Couette with finite aspect 
ratio, i.e., if the cylinders are covered with rigidly ro- 
tating end-plates (insulating or conducting), the Ekman- 
Hartmann layer also develops. Naturally influence of the 
vertical walls introduces additional important effects and 
direct quantitative comparison with the previous work is 
not possible. We must take into account that the fluid 
which was ejected due to the Ekman blowing mechanism 
must eventually get back due to the conservation of mass 
and finiteness of the container. Nevertheless we will show 
that the rotating end-plates induce the Hartmann current 
which can change the global properties of the flow. 

Let us introduce the parameter a which measures the 
overall importance of the magnetic field, 

d E Ha 
a = —= — = , (6) 

V2d H \/2ReM 

where c?h = -DHa -1 is the Hartmann depth. For the Ek- 
man depth, as a measure of the uniform rotation we use 
O out . The magnetic effects start to be significant when 
a > 1, in the limit a — > we have the classical Ekman 
layer and for a — > oo the classical Hartmann layer. We 
notice that for slow rotation corresponding to Re of order 
O(10 2 - 10 3 ) and Ha of order O(10), a « 1, and therefore 
we expect the magnetic fields to be important for many 
laboratory experiments. 

A. Insulating end-plates 

First we consider a case when both the cylinders ro- 
tate with the same angular velocity f2i n = fi ut = 100, 
i.e., p, = 1.0, and the rotational profile fl} is flat. The 
aspect ratio is T = 10 and the insulating plates rotate 
with angular velocity slightly different than the cylinders, 

fiend = 90. 

Figure [1] shows how the axial velocity u z and the axial 
current j z change with distance z from the plates, for 
different strength of the applied magnetic field. It can 
be seen that the axial velocity and the axial current de- 
crease for stronger magnetic field. The explanation is as 
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FIG. 1: Structure of the Ekman-Hartmann layer for the MHD Taylor-Couette flow with Q- m = Q ou t = 100, F = 10, insulting 
end-plates rotate with Q cn d = 90, z is the distance from one of the plates. Different values of the magnetic interaction 
parameters correspond to: Ha -» 0.0 (a -» 0.0), Ha = 3.0 (a = 0.2), Ha = 10.0 (a = 0.7), Ha = 30.0 (a = 2.1), d E = 0.0173 is 
the Ekman layer thickness. 



follows. The vertical shears in ur and u z produce cur- 
rents which together with axial field generate body forces 
acting against the shears. Since the radial flow must van- 
ish at the boundaries as well as it vanishes far away from 
them, the effect is to reduce the ur and, due to mass 
conservation, u z . Therefore the external axial magnetic 
field inhibits the Ekman blowing (which is completely 
suppressed when a-*oo) and makes the boundary layer 
thinner. The azimuthal flow uj,, on the other hand, is 
forced to have different values at the boundaries and far 
away from them, thus the shear can be decreased only in 
the region close to the boundary. 

These results are in a good agreement with the linear 
solution [1H for the case of the infinite, rotating plate 
and Pm — > (the agreement for other quantities like ur, 
i*0 is pleasing as well). We notice that in the radially 
unbounded case. u z and j z are independent of R which 
can not be true for the enclosed Taylor-Couette system. 
The values presented in Fig. [1] are computed for R = 
Ri n + D/2, in the middle of the gap, so that the influence 
of the rotating cylinders is smallest. 

We point out that the induced axial current j z (the 
Hartmann current) exists outside the boundary layer. 
This is not the case for non-rotating Hartmann bound- 
aries. For unbounded flow this current quickly converges 
to an asymptotic constant value, but for the case of flow 
between two plates, or for the enclosed cylinders, it can 
not be true and currents induced by both end-plates must 
eventually interact. When we consider a system symmet- 
ric in the z direction, i.e., when the two plates rotate in 
the same manner, the induced j z have the same strength 
but opposite signs and they eventually meet turning into 
the radial direction (and consequently j z = in the mid- 
dle of the container for the symmetric boundary condi- 



tions). 

We have varied < f2 e nd < ^out for constant Qi n = 
^out = 100, similarly we considered < £l m = f2 ut < 
100 for £l C nd = 100 to get values of the Ekman/Hartmann 
blowing and suction when the difference between cylinder 
and end-plates rotation is large. The agreement with pre- 
vious nonlinear calculation for the infinite plate is quite 
good, j25[. The dependence of the induced mass flux and 
the current on the strength of the magnetic field as well 
as on the relative fluid/end-plate rotation has the same 
character. 

If the flow is vertically bounded by two plates, as 
for the Taylor-Couette system, three essentially differ- 
ent regions can be distinguished: the Ekman-Hartmann 
layer, a magnetic diffusion region and a current-free re- 
gion ( [111 H3] ) • I n the magnetic diffusion region (MDR) 
the axial Hartmann current must be reduced to zero be- 
fore it reaches the current-free region and, by continuity, 
it is turned into radial direction. This radial perturbation 
current interacts with axial magnetic field and results in 
accelerating (for negative jr and positive B z ) or deceler- 
ating (for positive ]r) electromagnetic body force. 

The MDR arises since the Ekman-Hartmann layer it- 
self is incapable to force the current to satisfy the exterior 
boundary conditions. It constantly grows in time and it 
quickly dominates the whole space between the plates. 
Moreover, when considering the small Pm limit, the 
MDR instantly becomes spatially uniform and infinitely 
thick even for one bounding plane and the current-free 
region does not exists at all, [27] . 

Consequently in our enclosed MHD Taylor-Couette 
system with Pm — > we have relatively thin Ekman- 
Hartmann layer close to the plates, whereas the fluid in 
the major part of the container forms the MDR in which 
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the axial Hartmann current changes into radial one. We 
underline here that this is true for perfectly conducting 
walls, since such radial boundary conditions assure us 
that the current can penetrate the cylinders. The sit- 
uation would be rather different with insulating radial 
boundaries. 



B. Conducting end-plates 

For highly conducting plates the induced current 
drawn into/from the plates is much stronger than the 
current induced in the layer for insulating boundaries. 
The Ekman-Hartmann layer itself is nearly unaffected 
by conductivity of the plates as are the velocities and the 
currents within this layer. However, due to constant mag- 
netic field perturbation there exists an additional electric 
current of order 2$ which is induced by the conducting 
boundaries, <f> = eyOplaW^) an d e characterizes the rel- 
ative conductance of the fluid and the thin plates, [lj|. 
Moreover, in the MDR this current increases fluid veloc- 
ity by a factor of 2a 2 $. 

Figure [2] shows how the radial current jn in the middle 
of the container changes with conductivity of the end- 
plates. The difference between the perfect insulator and 
the perfect conductor is almost one order of magnitude 
even for so slow rotation. We find that for an MR1 ex- 
periment it is crucial to use insulating plates in order to 
minimize this undesirable current. 



IV. THE INFLUENCE OF THE HARTMANN 
CURRENT 



We notice that the force due to the radial current and 
the axial magnetic field, Ha 2 jfle/? x Bq/Bq, enters the 



momentum equation (|5a|) for component. Formally 
the force is equivalent to applying an azimuthal pressure 
gradient d^p ^ 0. A flow between rotating cylinders 
with non-zero d^p is usually referred as the Taylor-Dean 
flow, [28|. Its rotational profile £!d is a superposition of 
the circular Couette profile {T]) and the steady flow 
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with 



n + e(c + d/R 2 + lnR) 



R? n ln(i? in ) — R 2 ut ln(i? out ) 



E>2 



Rl 



d = 



^out -^in 

— d<pp. 
pv 



(7) 

(8) 

(9) 
(10) 



The pressure gradient can be realized by an external 
pumping mechanism or, like in the discussed case, by 
the Lorentz force resulting from the induced current and 
the axial magnetic field. 

Let us introduce a parameter (3 describing Taylor-Dean 
flows, the ratio of average pumping velocity to the rota- 
tion velocity 



insulating 




6V„, 



^in^in 

where Vm is the average pumping velocity 



(ii) 



i 



Ho 



V m = - I [e (c + d/R 2 + InJJ)] dR 

L> JR in 

Rout (1-V 2 ) 2 -477 2 (lnr)) 2 



2pv 4(1 - ?7)(1 - r) 2 ) 



(12) 



[29(. The basic question arises whether the resulting 
pumping due to the radial current and the axial field 
can bring the flow into an unstable regime. 



A. Hartmann current generated by the end-plates 



FIG. 2: Radial current j R (R = R in /D + 1/2, z = T/2) in 
the middle of the gap for MHD Taylor-Couette flow with 
Ha = 10, fiout = n in = 200, Sl cnd = 202, V = 10. The up- 
per line represents insulating plates, the bottom line perfectly 
conducting ones, and in between is the intermediate case for 
different values of the relative conductance. 



The structure of the Ekman-Hartmann layer changes 
with parameters such as rotation rates or strength of the 
magnetic field. Here, however, we will concentrate on the 
flow in the bulk of the container so that only currents 
and velocities which leave the layer are important. We 
analyze hydrodynamically stable flow with jl = 0.27 at 
the aspect ratio T = 10 with rigidly rotating end-plates. 
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1. End-plates rotating with f2 ut and fii n 

First we consider cylinders covered with rigid, per- 
fectly conducting plates rotating with the angular ve- 
locity equal to that of the outer cylinder fi cn d = fi ou t- 
We choose conducting lids so that the induced current 
is much stronger and its influence on the flow is more 
evident. 



A B C D 




Re=100 Re=200 Re=650 Re=850 



FIG. 3: Contour lines of stream function for different flow pa- 
rameters, the left edge of each panel denotes the inner cylin- 
der, the right edge the outer one, solid lines correspond to 
clockwise fluid rotation. Fhe end-plates are attached to the 
outer cylinder, tt cnd = fi out , ft = 0.27, T = 10. "A", "B" are 
for conducing plates and Ha = 3, "C", "D" for insulating ones 
and Ha = 10. 

When the plates rotate with fi cn d = fiout, the Ekman 
circulation is clockwise and the corresponding Hartman 
current has the positive sign, i.e., close to the inner cylin- 
der it leaves the Ekman-Hartmann layer with j z > 0, con- 
sequently the radial current also has positive sign. Fig- 
ure [3] A-B displays a flow with conducting plates and a 
weak axial magnetic field applied, Ha = 3, for two differ- 
ent Reynolds numbers. The rotation ratio is fx = 0.27, so 
that the Couette flow is hydrodynamically stable, how- 
ever we notice that when Re is large enough the flow 
changes significantly and the Taylor vortices can be ob- 
served. 

This phenomenon can be explained as follows: for a 
constant Ha, increasing of the rotation rate leads to the 
stronger Hartmann current drawn into the flow, therefore 
the corresponding pumping (3 due to jne^ x Bo increases 
and for certain Re it reaches a critical value f3 c , so that 
the instability develops. 



If the perfectly conducting ends are replaced with in- 
sulating ones the induced current is much weaker. When 
the imposed magnetic field has strength such that Ha = 
3, the pumping is too small to make the flow unstable, 
regardless of the Reynolds number. However, when the 
magnetic field is stronger, Ha =10, for sufficiently high 
rotation rates the vortices can also be seen, Fig.[3]C-D. 

It is known that stronger axial magnetic field has a 
stabilizing effect even on a hydrodynamically unstable 
flow [e.g. [3fJ| . Besides that, the Hartmann current in- 
creases with the amplitude of the magnetic field only until 
a certain point is reached. When the magnetic interac- 
tion parameter becomes a « 2.5, increasing Ha does not 
further increase the Hartmann current, [24j . For these 
reasons it is clear that when the imposed magnetic field 
is strong enough the instability described above will not 
occur. Indeed, it has been checked that for conducting 
plates, Re = 200 and the magnetic field with Ha = 20 
there are no Taylor vortices, although the rotational pro- 
file is significantly changed when compared to the non- 
magnetic situation. 

If rigidly conducting end-plates are attached to the in- 
ner cylinder, so that fiend = fiin> the Ekman circulation is 
counter-clockwise (Ekman suction) and the correspond- 
ing Hartmann current has a negative sign, so that the 
parameter [3 is positive. From Fig. [4] E-F we see that, 
analogously to the case "C"-"D", if the rotation is suffi- 
ciently fast the resulting (3 reaches critical value and the 
flow becomes dominated by the vortices. 

Similarly when the insulating plates are used, the axial 
magnetic field with Ha = 3 is too weak to generate suffi- 
ciently large fj. When stronger field is applied, Ha = 10 
it is possible to observe the instability ("G"-"H"). 



2. The rotational profile 

As mentioned above, if the induced radial current has 
the same sign as the axial magnetic field, the azimuthal 
velocity of the fluid is decelerated, if the signs are oppo- 
site the flow is accelerated. The discussed instability is a 
centrifugal one and is simply due to change in rotational 
profile of the fluid. Let us use a Rayleigh discriminant for 
stability, C = d R {R 2 n)/(Rfl), the flow is stable if ( > 0. 
Fig. O shows the radial dependence of £ in the middle of 
the gap (z = T/2) for the two cases labeled as "B" and 
"F" . 

We notice that the vortices concentrate in region where 
£ is negative, i.e., where the Rayleigh criterion is not 
fulfilled. This instability has essentially local character, 
so it is not possible to define any specific critical Reynolds 
number whose crossing would lead to some exponential 
grow in the whole container. For the conducting plates 
and fiend = fiout, there exists Re between 100 ("A") and 
200 ("B") for which only a part of the container would 
be filled with the vortices. 
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Re=100 Re=150 Re=700 Re=1200 



FIG. 4: The end-plates are attached to the inner cylinder, 
Q, cn d = f2i n . "E", "F" are for conducing plates and Ha = 3, 
"G", "H" for insulating ones and Ha = 10. Note that for the 
insulating endplates Re is order of magnitude larger. 



B. Linear stability of current-induced MHD 
Taylor-Dean flow 

In order to predict the onset of the instability discussed 
above we analyze the global stability of MHD Taylor- 
Dean now for our parameters. For the nonlinear simula- 
tions we can estimate the pumping due to the azimuthal 
pressure gradient just by setting (Vp)^ to Ha 2 j_R, see 
Eq. (pja|) . Generally (Vp)^ and jr change with radius 
like R" 1 . However, due to the presence of the plates, 
for jn this is true only far from the vertical boundaries 
and here the value of jr is taken at R = R ln / D, z = T/2 
(note that for our perfectly conducting boundaries the 
current penetrates the cylinders and for a steady state 
it is largest at R = R ln /D). In this way we obtain the 
parameter (3 associated with the enclosed MHD Taylor- 
Couette for the given boundary conditions "A"-"H" , and 
then it can be compared with a critical value f3 c obtained 
from the linear stability analysis. 

Consider now the axisymmetric MHD Taylor-Dean 
flow for infinitely long cylinders governed by the Eqs. 
It admits the basic solution = R£Id with ur — u z — 
= bj, = and the imposed axial magnetic field Bq. 
The perturbed state is u' R , RQn + u'±, u' z ,b' R , b'i, B + b z . 

After developing disturbances into normal modes we 
seek solutions of the linearized MHD equations in the 
form similar like in [3fj| . [3lT | . An appropriate set of ten 
boundary conditions is needed in order to solve the sys- 
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FIG. 5: t(R,z = T/2) calculated for MHD Taylor-Couette 
flow with conducting end-plates attached to the outer cylinder 
("B", Re = 200) and the inner one ("F", Re = 150). The 
Rayleigh criterion yields £ > for stability. 

tern, these are the no-slip boundary conditions for the 
velocity u' R — u'^ — u' z — and perfectly conducting for 
the magnetic field dub'^ + b'^/R = b' R = at the both 
cylinders. We will only consider stationary marginally 
stable modes. 

The homogeneous set of equations together with the 
boundary conditions for the walls determine an eigen- 
value problem of the form L(/t, t), k, m, Pm, Re, Ha, 0) — 
0. The variables are approximated with finite difference 
method on a grid typically with 200 points. The numer- 
ical code used to solve the problem is identical to that 
used in [H 

For the current axisymmetric (see, however, [34j ) study 
we set parameters m — 0,fj — 0.5, fx — 0.27, Pm = 10~ 6 , 
then for given Ha and Re we look for minimal value of 
\j3\ leading to the instability (the value for which the 
determinant L is zero). Since (3 is directly proportional 
to the azimuthal pressure gradient, and therefore to the 
radial current, the resulting critical (3 C determines the 
minimum value of the radial current for which the Taylor- 
Dean flow becomes unstable. 

Figure [S] shows marginal stability lines for the MHD 
Taylor-Dean flow for different values of the imposed axial 
magnetic field, for both positive and negative values of 
(3. We notice that much larger values of \(3\ are needed 
for stronger axial magnetic fields since the field plays a 
stabilizing role. 

The labels "A"-"H" refer to MHD Taylor-Coutte flows 
presented in the previous section. E.g., "A" refers to the 
flow with Re = 100, Ha = 3, fi = 0.27 with perfectly con- 
ducting end-plates attached to the outer cylinder. The 
induced current jn is such that the corresponding f3 due 
to the Lorentz force denotes a stable flow. If the Reynolds 
number is increased, the critical value (3 C (for Ha = 3) is 
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FIG. 6: Critical values of the pumping/rotation ratio j3 c for MHD Taylor-Dean flow with ji = 0.27, for different Reynolds 
number Re and strength of the magnetic field Ha. The case for negative (positvie) f3 corresponds to the azimuthal pressure 
gradient due to positive (negative) radial currents interacting with the axial field. On the left panel the unstabel region lies 
below the line, on the right one above the lines. The letters "A"-"H" represent states displayed in Figs. [3] [4] 



reached and the instability develops - label "B" . 

V. SUMMARY 

Gilman and Benton [24[ have shown with a linear the- 
ory that in vicinity of a rotating plane which serves as a 
border for rotating conducing fluid there develops the 
Ekman-Hartmann layer if Opiate 7^ Ofluid and an ax- 
ial magnetic field is applied. The most important fea- 
ture of Ekman-Hartmann layers is their ability to in- 
duce both mass fluxes and electric currents in the re- 
gion outside the boundary layer. If Opiate < Ofl u id these 
fluxes are directed outwards the layer ( "blowing" ) ; when 
Opiate > Ofl u id towards the layer ("suction"). For the 
conducting plates the fluxes are much stronger since ad- 
ditional currents are drawn from/into the plates. 

Outside the Ekman-Hartmann layer exists the mag- 
netic diffusion region, in which the electric current has 
only radial components. The current, together with the 
axial magnetic field, produces an electromagnetic body 
force acting on the fluid. 

We have shown in this paper that similar effects arise 
for the MHD Taylor-Couette flow when the rotating 
cylinders are bounded by two rigidly rotating end-plates. 
Near the plates the Ekman-Hartmann layer forms and, 
consequently, there exists the Hartmann current which 
penetrates bulk of the fluid. In the presence of an ax- 



ial magnetic field such problem can be compared with 
the Taylor-Dean flow - a flow between, possibly rotating, 
cylinders which is additionally driven by an azimuthal 
pressure gradient. 

We find that under certain conditions the resulting flow 
becomes unstable, Taylor vortices can be observed and 
the rotational profile is significantly different from the 
standard Couette solution Oo. The instability has essen- 
tially a centrifugal character as the Rayleigh criterion is 
locally vialoted. This is an undesirable effect from the 
point of view of an MRI experiment. In such experiment 
it is necessary to obtain a state resembling Oo in the ma- 
jor part of the container, for parameters characterizing 
stable MHD flows. It is necessary to take into account the 
magnetic effects induced by the plates so that the MRI 
can be clearly identified rather then any other instability. 

The fluxes induced in the Ekman-Hartmann layer are 
a direct consequence of a shear close to the boundaries. 
Exemplary methods of reducing the shear have been pro- 
posed in For rotation rates characterized by Re of 
order O(10 3 ) all the effects can by significantly reduced 
by allowing the end-plates to rotate independently of the 
cylinders [32j ] . Since for O on d = Oj n there is the Ekman 
suction, and for O on d = O out the Ekman blowing, there 
exists O out < O en d < 0; n for which the generated mass 
and charge fluxes are minimal. Alternatively, one can 
divide the plates into independently rotating rings [33j |. 
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